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Abstract 



We obtain the exact solution of a general Hubbard chain with kinetic 
energy t, bond-charge interaction X and on-site interaction U with the only 
restriction t = X. At zero temperature and half filling, the model exhibits 
a Mott transition at U = 4t. In the metallic phase and near half filling, 
superconducting states are part of the degenerate ground state and are favored 
for small U if the system is slightly perturbed. 
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The exact solutions, particularly those obtained using the Bethe ansatz, have brought a 
very important progress in the understanding of strongly correlated systems. However the 
conditions for integrability using the Bethe ansatz are very restrictive and only a limited 
class of realistic models can be solved with this technique [Q. Due to the importance of 
the exact solutions in clarifying the effect of different physical ingredients and as a test of 
approximations, the search of exact solutions has been recently extended to other models 
and techniques, in spite of the fact that in some cases the model or the parameters are rather 
unrealistic 

The model we consider is a particular case of the following Hamiltonian: 

H = Hu + Ht = UY, rii^^ii + 

i 

<ij>(T 

- rij„) + (1 - ni„)nj„] + tBBniaUj^} . (1) 



H has been derived as an effective one-band Hamiltonian for the description of cuprate 
superconductors Similar models including in some cases the nearest-neighbor repulsion 
V have been studied by several authors PJ5|J5|-p!3| . If Iaa + tBB — '^^ab = the three-body 



term of Ht vanishes, and H reduces to the model considered by Hirsch and Marsiglio, in the 
framework of their theory of "hole superconductivity" |jlO|- Following Ref. f^, we call the 
coefficients of the one-and two-body parts of Ht as Iaa = ~t and tAB^tAA = X respectively. 
In the weak-coupling case < X << t, a standard BCS-type mean-field approximation 



and a renormalization-group analysis in the one- dimensional (ID) continuum-limit theory 
[131 , show that a small positive X gives rise to an effective atractive interaction for a particle 



density n > 1, while this interaction is repulsive for n < 1, and vanishes at half filling. 
This situation cannot be extended to the case X = since for these parameters {Iab = 
tAA + tsB = Hf is symmetric under an electron-hole transformation and the physics for 
densities n and 2 — n should be the same. Thus, it is of interest to study this case. This is 
one of the goals of this Letter. Strack and Vollhardt studied the model for these parameters 
(including V) at half filling and argued that this case correspond to a physically relevant 



range of parameters 0. 

The study of the Mott transition also makes the case t^B = appeahng, because of the 
supression of antiferromagnetic correlations. This avoids the problem of having to distinguish 
between a Mott insulator in which the particles become localized as a consequence of strong 
on-site repulsion and an antiferromagnetic insulator, in which a weak interaction opens a 
gap in a nested Fermi surface. The latter is the case of the Hubbard model in bipartite 
lattices. Studies of the Mott transition in these cases are restricted to the paramagnetic 



phase |l^-|l6i- Other studies have taken nonbipartite lattices or systems in which the 



noninteracting Fermi surface has no nesting ||T8|,[T9|. In the large U limit, the model of Eq.(l) 



becomes equivalent to a generalized t — J model with hopping tAAitBs) for n < l{n > 1), 
correlated hopping t\^/U, and antiferromagnetic exchange interaction J = 4:t\^/U which 
vanishes for tAs = 0. 

In this Letter we obtain the exact solution of Hamiltonian (1) for a chain with open 
boundary conditions under the only restriction tAB =| ^aa | — | ^bb |= 0. We also discuss 
the effect of a finite tAB on the basis of our Lanczos results for finite chains. Strack and 
Vollhardt obtained the exact ground state for Ibb = —^aa = t, for arbitrary dimension 
including the nearest-neighbor repulsion V, but only for n = 1 and two regimes of parameters 
in which all particles are static in the ground state [§]. In ID and for V = we are able 
to obtain all eigenstates for arbitrary filling, particularly in a third regime of parameters 
in which the dynamical part of the Hamiltonian plays an important role in the ground 
state. 

The exact solution of the model is greatly facilitated by its symmetries. In any dimension 

for tAB = 0, [Ht, Hu] = and the number of doubly occupied sites is conserved 0. Also, as 

in the case of the model of Essler, Korepin and Schoutens 0, for tAB = 0, Ht commutes not 

only with the total spin, but also with the following generators of another SU(2) algebra: 

L L L 

V = Y1 c^i^iT, = Y1 4Av Vz = J2(2-Yl clcia), (2) 

i=l i=l i=l CF 

where L is the number of sites. This allows us to construct eigenstates of minimum energy 



which posess off-diagonal long-range order for sufficiently small values of U and | n — 1 | . 

The solution of the chain is obtained mapping Ht into a tight-binding model of spinless 
fermions. To obtain this mapping it is convenient to write if in a slave-boson representation. 
We represent the four possible states at site z :| >, | >, c||cj| | >, by e| | >, //^ | 
>,c/| I > (pictorially o, | or | and •) respectively, using two bosons to represent the 
empty (o) and doubly occupied (•) sites and two fermions (| and [) to describe the singly 
occupied sites. The Hamiltonian takes the form 

H = UJ2 4di + tAA Yl fjafiaeUj 
<ij>^ <ij> 

+h.c.), (3) 

with the constraints e|ej + dldi + J2a fLfi^r = 1- When t^B = 0, the numbers = 
J2i fLfi(T: = Hi^l^i and Nd = J^idldi are separately conserved. Note also that in a 
bipartite lattice, changing the phase of the bosons Cj or di by (-1) in one sublattice changes 
the sign of t^A or tss respectively. Thus we can choose these signs arbitrarily. Taking 
—tAA = tBB = t > as in Ref. 0, Ht takes the form: 

Ht = -t E [flfUeUj + did,) + h.c] (4) 

<ij>(T 

In a chain with open boundary conditions also the order of the bosons and that of the 
fermions along the chain are separately conserved: Ht permutes the order of a fermion and 
a boson which are nearest neighbors, but two bosons or two fermions cannot be permuted. 
For a given number of fermions Nf = N-^ + Ni, let us numerate the L sites, Nf fermions and 
Nf, = Ne + N(i = L — Nf bosons with similar sequence (for example from left to right) using 
the labels i,j and m respectively. Then, any state with definite number of particles on each 
site can be written as: 

I >= n [B{m)el^^^ + (1 - 5(m))4^)] 

m=l 

Nf 

xnTO)4)T+(i-^o'))4)jio>- (5) 
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Here i(m) is the position of the mth boson in the sequence (its inverse, defined on the 
set of sites for which rife, = ejcj + = 1 is simply m{i) = Y4=inbi), and has a 
similar meaning for the fermions. B{rn) = 1 if the mth boson is an "empty" one and 
zero otherwise. Similarly in terms of the spin of the fermions F{j) — 1/2 + 'S'?^^). The 
products are ordered throughout with increasing labels to the right. As an example the 
state I il^i >= o IJ, o« IIJ, • ... and any other state | ipi^ > such that < ipi \ Hf \ ipi^ >^ 
have B(l) = B(2) = 1, B(3) = S(4) = 0, F(l) = F(3) = F(4) = 1 and F(2) = F(5) = 0. 

Due to the properties of Eq.(4) and the open boundary conditions, the ID model has 
an extremely rich symmetry structure, including L SU(2) symmetries which are the local 
versions of those previously mentioned. There is one usual spin SU(2) algebra related to 
each of the Nf fermions and a "local pairing" SU(2) algebra related with each boson. As an 
example it can be easily verified that {Htel^^-^di(m)-4(m)di(^)^t) \ ipi >= 0, where el^^^d^m) 
is a raising operator. Thus one can separately diagonalize Ht in each subspace of definite 
values of B(m) and F(j). For fixed Nf there are 2^ subspaces and the size of each one 
is (Nf)- The raising and lowering operators establish a one to one correspondence between 
each state of one of these subspaces and the corresponding one of another subspace and Ht 
takes the same form in all these subspaces. In the subspace of highest weight of all SU(2) 
algebras (all i?(m) = F{m) = 1), the solution of Ht for given Nf is easily obtained. The 
eigenstates, written in the original representation have the form: 

I € >= n 4,T 1 >, 4, = {Y—y/'T.Mki)ci (6) 

i=i ^ i 

where the possible values of k{L + l)/'!r are positive integers. These eigenstates can be 
extended to any values of B{m) and F{j) using the lowering operators: 

I V'e >= n {rifi[FUi) + (1 - FUi))4iCir]+ 

1=1 

+(1 - nf,)[B{m,) + (1 - B(mO)c,Vlj} | >, (7) 
where Ufi = nj(2 - n^), rii = J2a cla^iaji = El=i rifi and rui^i- ji. 
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Eqs. (6) and (7) also describe all the eigenstates of H = Ht + Hu- The latter term 
reduces the degeneracy to 2^-' (j^^) and adds UNa to the energy. 

For each particle density n, the ground state of H is obtained minimizing the energy 
with respect to the density of doubly occupied sites d = Nd/L and taking the lowest Nf 
values of k in Eq.(6), with the constraint nL = 2dL + Nf. The result is very simple. In 
the thermodynamic limit three regimes can be distinguished depending on the values of U/ 1 
and the particle density n. Also three regions of values of U/t can be separated (For n = 1 
and \ U \> 8t the ground state was obtained previously by Strack and VoUhardt |^): 
a) U > At. This region lies inside what we call regime I: for n < 1 the physics is the same 
as that of a spinless model. The ground state expectation value < Hu >= and: 

2t 

d = 0, e{n) = sin(n7r), (8) 

TT 

where e(n) is the energy density. For n > 1, from electron- hole symmetry d = n — 1, e{n) = 
U{n — 1) + e(2 — n). For n = 1, < Ht >=< Hu >= and the system is an insulator with 
energy gap U — At. 

h) U < —At. This region coincides with regime II. Here (for an even number of particles) 
all particles are paired, all pairs are static (< if^ >= 0) and: 

d = n/2, e{n) = Un/2. (9) 

c) —At < U < At. In this region there are two critical densities rii and n2 defined by: 
rii = (I/tt) arccos(— t//4t) and ni < 1 < ^2 = 2 — rii. For n < ni or n > n2 the physics 
corresponds to regime I and the ground state and its energy was described above. Instead, 
for Til < n < n2 the system is inside regime III. This regime is the only one in which empty, 
single and double occupancy at any site is possible, and the competition between Ht and 
Hjj is apparent in the ground state. The double occupancy and energy are given by: 

d = e(n) = Ud- —{IQe - U^Y'^ (10) 

2 27r 

In regimes II and III the system is at the borderline of phase separation and also of 
superconductivity. Eigenstates with off-diagonal long-range order (ODLRO) are part of 
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the degenerate ground state. To show this, let us take an eigenstate \ ipg > of the form 
of Eq.(7), with A^^^ doubly occupied sites, which belongs to the ground state. The state 
I ip >= 7]^"^ I ipg > with 1] given by Eq.(2), is clearly different from zero (it is obtained 
from I tpg > putting all B{m) = 1 in Eq.(7)) and is also an eigenstate of Ht with the same 
eigenvalue as that of | ipg >. Also | > is a highest-weight state of the ?7-pairing SU(2) 
algebra (Eq.(2)). Similarly the state | ipN^ >= (v^)'^'^ \ ip > is an eigenstate of Ht with the 
same eigenvalue, and an eigenstate of H with the same energy as the original state \ ipg >■ 
In Ref. 1^, it is shown that | ipN^ > the thermodynamic limit (L — oo with d = Nd/L 
constant) has ODLRO if 7^ and l + d-n = N^/L ^ 0. 

The model has a metal-insulator transition at Uc = 4t. The four-boson theory of Kotliar 
and Ruchenstein |]15[ in the mean-field approximation gives Uc = IGt/vr in good agreement 



with the exact value. The approximation also gives a reasonably accurate Uc for the infinite- 
dimensional Hubbard model [|1^ . 



The form of the Hamiltonian in the representation of Eq.(3) suggests that addition of a 
small t^B such that it can be treated in second-order perturbation theory, introduces anti- 
ferrogmagnetic correlations between nearest-neighbor fermions and allows the permutation 
of nearest-neighbor bosons d and e, increasing their mobility and favoring superconductiv- 
ity. We have solved numerically the model for tsB = —tAA = ^,tAB = 0.2andL = 10. For 
1/2 < n < 1, the model exhibits phase separation for U > Ug with f/^ ~ 1 for n ~ 3/4 
and f/5 = for n = 1, while for U < Ug the system behaves as a Tomonaga-Luttinger 
liquid (TLL) |21[]. For n < 1/2 the TLL behavior is observed for all values of U. Within 
the TLL regime, the evaluation of the compressibility, the Drude weight and the spin and 



charge velocities allowed us to derive the correlation exponent Kp ||2T| . The resulting values 
indicate that the dominant correlations are the superconducting ones for 1/2 < n < 1 and 
the charge-charge ones for n < 1/2. 

In this Letter we have solved exactly a Hubbard chain including bond-charge repulsion 
for a particular value of the latter. The model displays a Mott transition at half filling and 
in two regimes of parameters the ground state contains superconducting states. Numerical 



results show that superconductivity is favored by a small perturbation for not too large 
on-site Coulomb repulsion. 

One of us (L.A.) is supported by the Consejo Nacional de Investigaciones Cientificas y 
Tecnicas (CONICET), Argentina. A. A. A. is partially supported by CONICET. 
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